## 6.1. Crie seus dados ##
#crie dois conjuntos de dados
A<-rnorm(10,mean=6,sd=3)
A
B<-rnorm(10,mean=7.5,sd=3.2)
B
#anlises exploratrias
summary(A)
summary(B)
hist(A)
hist(B)

#diferena das mdias
diferenca<-data.frame(variaveis=c(A,B),fator=factor(rep(c(1,2),each=10)))
diferenca
boxplot(diferenca$variaveis~diferenca$fator,outline=FALSE)
mean(A)-mean(B)
mean(B)-mean(A)
dif<-abs(mean(A)-mean(B))
dif
hist(diferenca$variaveis,freq=FALSE,xlim=c(0,14))
curve(exp=dnorm(x,mean=mean(diferenca$variaveis),sd=sd(diferenca$variaveis)),from=0,to=14,col="red",add=T)

##simulao
source("simula.r")
X11()
smla.diferenca<-simula(A,B,1000,"bi")
# as mdias no so diferentes

smla.maior<-simula(A,B,1000,"uni")
# As medias do B nao sao maiores do que do conjunto A

t.dif<-t.test(A,B,"two.sided")
t.dif
t.maior<-t.test(A,B,"greater")
t.maior
##Os resultados corroboram com os resultados das simulaes.